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Abstract. We generalize the Maximum Likelihood-type method used to 
study cross correlations between a catalog of candidate astrophysical sources 
and Ultrahigh Energy Cosmic Rays (UHECRs), to allow for differing source 
luminosities. The new method is applicable to any sparse data set such as 
UHE gamma rays or astrophysical neutrinos. Performance of the original and 
generalized techniques is evaluated in simulations of various scenarios. Applying 
the new technique to data, we find an excess correlation of about 9 events between 
HiRes UHECRs and known BLLacs, with a 6 X 10~^ probability of such a 
correlation arising by chance. 



1. Introduction 

Correlation studies play a fundamental role in establishing or ruling out candidate 
sources for rare events such as the highest energy cosmic rays and UHE gamma 
rays. When the number of events is low, it is necessary to use sensitive methods 
to be able to identify the sources and quantitatively assess the possibility of an 
incorrect identification. The purpose of the present work is to evaluate and improve 
these methods. If the hypothetical sources are commonplace rather than rare, so 
that the typical separation between candidate sources is not large compared to the 
directional uncertainty of individual events, then an entirely different approach must 
be used, requiring much larger datasets. Fortunately, the GZK horizon allows the 
surface density of source candidates to be reduced by increasin g the energy thr e shold ; 
this played an essential role in the recent Auger discovery ( Abraham et al.l . 120071 ) 
of a correlation between UHECRs and galaxies in the Veron-Cetty Veron Catalog 
of Quasars and AGNs. Furthermore, it is plausible that only unusual astrophysical 
objects can produce very rare events, in which case the candidate source catalog may 
naturally be sparse. 

Hints of excess correlations we re reported earlier between ultrahigh energy 



cosmi c r ays (UHECRs) a nd BL Lacs by lTinvakov and TkachevlpOOlLICorbunov et al.l 
12004"), 'Abbasi et al.' ('200^ ""-^ ^-^- — TTUT?r'T3„ „„^ „ „„„ „i,.„.^„.„ k„ 



Pierpaoli and Fa r rar (200 



and between UH ECRs and x-ray clus t ers b y 
The analyses in Tinvakov and TkachevI ( 20011 ) , 



Gorbunov et al.l (|2004l ). and IPierpaoli and Farraij (J2005r ) determine the number of 



correlated events within some given angular separation between UHECR and source, 
and calculate the "chance probability" of finding a correlation at the observed level by 
doing a large number of simulations with no true correlations. The ([Abraham et al. . 
20071 ) analysis in addition "scans" to find the optimal UHECR energy threshold and 



maximum source redshift. 
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In order to incorporate the experimental resolution on an event-by-event basis, 



Abbasi et alj (|2006h proposed a Maximum Likelihood-type procedure and applied it 



to studying correlations between UHECRs and BLLacs. This procedure (denoted 
the HiRes procedure, below) is motivated under the unphysical assumption that 
every candidate BLLac source has the same apparent luminosity. Even if BLLacs 
were standard candles with respect to UHECR emission, the BLLacs in the catalog 
have a large range of distances which would imply an even larger range of apparent 
luminosities, so one does not want to rely on such an assumption. Furthermore the 
validity of the method was not demonstrated via simulations. 

In the present paper we introduce a ML prescription which avoids the assumption 
of equal apparent source luminosity and allows the potential sources to be ranked 
according to the probability that they have emitted the correlated UHECRs. We test 
and compare both methods with simulations in a variety of situations. We find that 
the HiRes method gives the correct total number of correlated events even when the 
sources do not have equal apparent luminosities, as long as the numbers of events are 
sufficiently low and the candidate sources are not too dense or clustered themselves. 
We find that our new method performs better in these more challenging cases. In a 
final section, we apply the new procedure to BLLacs and x-ray clusters. 

2. Maximum Likelihood Approach for the Cross-Correlation Problem 

2.1. The HiRes Maximum Likelihood method 



In the HiRes Maximum Likelihood method of lAbbasi et al.l ( 20061 1 . the aim is to find. 



among N cosmic ray events, the number of events, n^, that are truly correlated with 
some sources, of which the total number is M . Hence, there are N — rig background 
events whose arrival directions are given by a probability density i?(x), which is simply 
the detector exposure to the sky as a function of angular position, x. For a true event 
with arrival direction s, the observed arrival direction is displaced from s according 
to a probability distribution Qi{x,s). Note that Q,: is governed only by the detector 



resolution. For the analysis given in I Abbasi et al.l ( 20061 ). Qi is taken to be a 2d 



symmetric Gaussian of width equal to the resolution ai, of the ith event. (Note that, 
throughout, the parameter a is related to ergs, the radius containing 68% of the cases, 
by ctqs = 1.51(7.) Dispersion due to random magnetic fields can be incorporated into 
the Qi by generalizing the variance to a^jf^^^.^^^ = aj^tector+'^^nagnetic: where amagnetic 
may vary with direction and event energy. In practice the magnetic dispersion is not 
well known and has to be treated as an unknown or argued to be smaller than a detector ■ 
The probability density of observing the ith event in direction x^ is 

" ■ Etii Q(xi, Sj)i?(Sj) N - n. 



and the likelihood for a set of N events is defined to be C{ns) = ni=i Pi{'^i)i which 
is maximized when Ug is the true number of correlated events. Since £ is a very 
small number, which depends on the number of events, it is more useful to divide 
C by the likelihood of the null hypothesis, i.e., Us = 0, to form the likelihood ratio 
TZ{ns) — C(ns)/C{0). The logarithm of this ratio is then maximized to obtain the 
number of correlated events, Ug. The significance of the correlation is evaluated by 
measuring the fraction JF of simulated isotropic data sets with as large or larger value 
of ln7^. 
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2.2. Extending method to differing luminosities 

If the source candidates for UHECRs were standard candles with known distances, 
or if the relative fluxes from different sources were known, we could generalize the 
above method by simply attaching the appropriate relative weight for each source and 
maximizing the likelihood ratio to obtain the number of correlations, n^. However, 
for the applications we have in mind we do not know the relative luminosities of the 
putative sources, nor do we expect them to be standard candles. Thus, to generalize 
the HiRes method to allow for sources with differing luminosities we assign an a apriori 
unknown number of cosmic rays, Uj, separately for each source, with ritot — X^^i ^j- 
The probability density generalizes to 

^M ^/ N n/ \ AT V^Af 



E,=i nj Q{x^,Sj)R{sj) iV^-Ej-1 ni 
NRs N 



P,(x,) = -'=' -;;-- " -' + — ^i^L^i?(x,), (2) 



where Rs = X]i=i ^i^j)/-^- ^^ above, we divide the likelihood C by the likelihood of 
the null hypothesis, C{{nj} ~ 0), to obtain the likelihood ratio 



N 



n^ll 



i=l 



T,f=injQ{^z,Sj)R{sj) i£. 



A^^o 



R{^z)Rs 



(3) 



Maximizing In 7^, with respect to the set {uj} of M numbers, determines the most 
probable values of the set {uj} of M numbers as well as the set {{InTVji} of A^ 
numbers. The numbers {{hiTZ)i} are figures of merit, providing information about 
how strongly correlated the individual cosmic ray events are to the catalog of sources, 
and allow us to rank the cosmic rays in order of their likelihood of being correlated to 
the source data set (similar information can be obtained in the HiRes approach from 
the contributions of individual cosmic rays to the sum in ln7^(ns)). 

A crucial difference between the generalized method and the HiRes method is 
that in the new method ritot = Yl ^i provides an estimate of all correlations, i.e., both 
true and random correlations, whereas the HiRes method yields only an estimate of 
the number of true correlations. In the generalized method one obtains Uj > if any 
cosmic ray event is close enough to any source, regardless of the degree of correlation 
between the two data sets elsewhere. For the HiRes method, the single optimized 
parameter, Us, will be greater than zero only if the degree of correlation between the 
two data sets are greater than what is expected for a random sample of cosmic rays 
(weighted by the detector exposure). 

We now turn to the estimation of Us with the new method. A crude estimate 
of the number of true correlations for the new method is ritot — nrand, where firand 
is the average number of correlations obtained when cosmic rays are uncorrelated to 
the data set of potential sources, i.e., cosmic rays drawn from an exposure- weighted 
isotropic distribution. 

A better measure of Us is summarized by the following equations, where the 
superscript in parentheses labels the "order" of refinement, giving successively better 
approximations to n^: 

M 



n 



^'^ = r\n^"^ - TirandiN)) (5) 



n 



f-\n'^"^-nrandiN-nW)), (6) 
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where Uj-andiN) is the average total number of correlations found with N events 
drawn from an isotropic distribution weighted by the detector exposure, and where / 
is the average fraction of true events that are recovered as correlated. This fraction 
will typically be slightly smaller than unity since under the assumption that true 
correlations are separated according to a Gaussian distribution, some events are 
separated too far from their sources to be accepted as correlated events. We measure 
/ by simulations. 

Up to this point, we have not specified whether {rij} are real numbers or integers. 
Indeed, the method allows either choice. A non-integer implementation has the virtue 
of encoding more information into the set of numbers, {rij}, while in the integer case it 
is necessary to also study {(lni?)j} in order to rank the sources in terms of correlation 
quality to the set of cosmic rays. An integer approach may be preferred as it provides 
a more concise answer as to which are the likely correlated sources. Moreover, in the 
case of strong correlations and ample statistics, the integer implementation is likely 
to offer a more intuitive result in terms of singlets, doublets, etc.. In the following 
sections we will use the integer implementation, unless noted otherwise. 

The procedure to maximize InTZ may depend on whether {rij} are non- integers 
or integers. In the present work, we use a commercial non-linear optimization package 
for the former case. One might expect to encounter problems when the number 
parameters to be fitted (M) exceeds the number of data points (N). However, in 
practice this never becomes an issue as long as Mna'^^^^ is small compared to the 
solid angle observed so the number of random correlations is not too large. Then, 
most potential sources are well separated from most UHECR events, and thus have 
rij = 0, thanks to the rapid cutoff of Q for s :^ a, where s is the separation between 
cosmic ray and source. For the integer case we first note that the exponential form of 
Q and the diluteness of the source data set will make it extremely rare for a correlated 
event to be ascribed to the wrong source. In other words, a cosmic ray event must 
be at almost the exact same (small) angular distance from two potential sources for 
there to be any ambiguity as to from which source the event originated. For the ith 
cosmic ray, this allows us to consider only the source j that maximizes the quantity 
(5(xi, Sj)i?(sj), and put that quantity to zero for all other sources. The maximization 
of equation [3] is then a matter of testing, for each source, whether \nTZ decreases 
when rij is increased from to 1. If it does, then rij — 0. If InTZ increases, then we 
repeatedly increase rij by unity, until InTZ decreases and wc find the correct number 
of correlations, rij, for that source. 

3. Simulation trials 

In vetting the two methods with simulations, we test their ability to correctly 
reproduce the number of true correlations on mock data sets. This allows us to 
explore the effect of large event and source densities, the effect of anisotropy in the 
source distribution, and the consequences of having incorrect event resolutions. 

3.1. Dilute source and UHECR data sets 

We begin by testing the ability to reproduce the number of true correlations in 
the simplest case of dilute, random sources. A first simulation is done using half 
the sky, uniform detector exposure -R(x), 156 randomly dist ributed sources and 271 



cosmic rays (the numbers relevant to the BLLac studies of iGorbunov et al.l (|2004l ) 
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Figure 1. The number of correlations vs. the detector resolution. Error bars 
(slighly separated horizontally for readability) contain 90% of the cases. Top: 
Random sources; Bottom: actual BLLac positions —30° < <5 < 90°. 



and lAbbasi et al.l ( 20061 )). Ten of the cosmic rays are Gaussianly aligned (a cosmic 
ray paired to a source, with angular separation according to the probability density 
Q, taken to be a 2d Gaussian of width cr), for various event resolutions. A second 
simulation uses the actual BLLacs source positions; these are more clustered than the 
random case. As shown in figure[Tl both methods reproduce well on average the correct 
number of correlations. The error bars shown include 90% of the 10k realizations, and 
are similar for both methods. As a increases, the variance in the number of found 
correlations increases rapidly. 



3.2. Dispersion in results 

To test the extent of agreement between the two methods in individual realizations, 
we generate 156 source positions at random on one hemisphere, align 10 cosmic rays 
to 10 of the sources and distribute 161 cosmic rays at random, then determine the 
number of correlations identified, with both procedures. Figure [5] shows the results of 
1000 repetitions, for cr — 0.4° and for 0.8°. From this exercise we conclude that the 
two ML methods do not generally agree for individual realizations, in spite of the fact 
that both methods reproduce the correct value in the mean. 
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Figure 2. Comparison of the extracted number of correlations for specific 
realizations using the two methods. 1000 realizations of ten cosmic rays aligned to 
sources. Left: For cosmic rays with event resolution a = 0.4° . The banding is due 
to using the integer implementation of the generalized method. Right: u = 0.8°. 



3.3. High source and UHECR densities 

When there are very large numbers of events and of potential sources, we expect 
the generalized ML method to perform worse than the HiRes method, beacause 
rv-'^> is obtained by taking the difference of two very large numbers, n'-"-' and ftrand- 
Moreover, as the source density becomes very large it becomes impossible to reliably 
distinguish the "contributing" sources. The total number of correlations becomes the 
only interesting quantity to calculate. Thus, only the HiRes method should be used 
for the case of very high event densities. However, the HiRes method also deteriorates 
at high densities, as shown in figure [3l 

5.^. Sensitivity to experimental resolution 

If the resolution of cosmic ray events are consistently over- or underestimated in a 
given data set, the extracted correlations will be incorrect. To test the sensitivity of 
the two methods to this problem we repeat the first type of simulations, but rescale the 
event resolution when aligning a cosmic ray to a source. Figure 2] shows the average 
number of correlations found, as a function of the amount by which a is rescaled, 
for the two methods. The new method is far less sensitive to incorrectly estimated 
resolution than is the HiRes method. 

3. 5. Clustering of sources 

As seen in the lower panel of figure [TJ spatial correlations within the data set of 
potential sources may skew the number of UHECR correlations found. In figure [5] 
we show the results of a simulation with clustering of potential sources. The figure 
shows the mean for 10k realization of 271 cosmic ray events with a = 0.4°, and two 
different scenarios for the correlation with source clustering, for 156 candidate sources. 
In both cases, candidate sources and CRs are distributed over one hemisphere; 100 
sources are placed in ten randomly positioned clusters with ten sources distributed 
around each cluster center according to a 2d Gaussian of width d, and the remaining 
56 candidate sources are placed at random in the hemisphere. In the first case, a 
randomly selected source in each cluster has one cosmic ray event Gaussianly aligned 
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Figure 3. Correlations found by the HiRes method with 10 true correlations, M 
sources and A^ events in 100 square degrees. The method's prediction becomes 
increasingly inaccurate, as the dilute approximation becomes less valid. 



to it and the remaining 261 cosmic rays are placed at random. In the second case, ten 
CRs are Gaussianly aligned with ten of the randomly placed source candidates not in 
any cluster and the remaining CRs are placed at random. As figure [5] demonstrates, 
the HiRes method significantly overestimates the true number of correlations if the 
sources are in clustered regions and underestimates it when the candidate sources 
show significant clustering but the UHECRs do not come from the clustered regions. 
By contrast the new method performs well in this test. 

4. Application to BLLacs and x-ray clusters 

In this section we apply both the HiRes method and the generalized method on two 
previously investigated data sets of suggested UHECR sources. 



4-.1. X-ray clusters 

The possibility that X-ray clusters are sources for UHECRs was investigated by 
Pierpaoli and Farran (J2005l ) using a binned analysis; a correlation w as found between 
AGA SA events and x-ray clusters taken from the NORAS catalog ( Bohringer et al.l . 
2000|). Here we study the correlation between X-ray clusters and two samples 
of UHECRs, both sep arately and combined. We use the 375 X-ray clusters 
(JBohringer et all . l2000f) in the region 0° < dec < 80°, and \b\ > 20°. For cosmic 
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Figure 4. Average number of found correlations as a function of the factor by 
which <j is reseated by when "Gaussianly aUgning" sources. 
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Figure 5. Sensitivity to clustering in source dataset, for UHECRs from (A) dense 
or (B) sparse regions. 



rays, we use the 143 HiRes events with E > 10"'^^ eV, and 36 AGASA events with 
E > 4:X 10^^ eV, in the same angular region. Note that in lPierpaoh and Farraij ( 20051 ) 
the authors restricted the correlation search to x-ray sources within the estimated GZK 
maximum distance for each UHECR event. Imposing the GZK r estriction reduces the 
numb er of found correlations in simulated random samples, but iPierpaoli and Farrar 
( 20051 ) found that it does not significantly decrease the number of correlations between 
UHECRs and x-ray clusters so that it increases the significance of the correlation. In 
this paper we have not taken the GZK distance restriction into account, as our focus 
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Table 1. Correlations between X-ray clusters and UHECRs. The integer 
implementation has been used for the generalized method. For comparison, using 
the non-integer implementation we find n'^^ = 3.3 and n(^) = 9.6 for the HiRes 
and AGASA event sets, respectively. 



Generalized ML Method HiRes Method 

UHECR N n(0) n(2) :f urr ^Fhb^ 

HiRes 143 22 5.2 0.3 1.5 0.3 

AGASA 36 24 9.9 610"^ 9.2 9 10"^ 

Both 179 37 12.7 0.06 4.6 0.1 



is on comparing the two ML methods. 

When applying the Maximum Likelihood methods to x-ray clusters, we need to 
take the angular size of the clusters into account, as they can be of the same order as 
the cosmic ray event resolution. Therefore we introduce an effective event resolution 
to be used in the calculation of Q, cTeff — vct^ + t^, where r is the angular radius 
of the cluster. When combining the data sets we calculate the detector exposure 
by i?(x) = X^i -^i-^i(^)/Si -^jj with i labelling the different data sets. For the 
generalized method we use the integer implementation. 

The results of the various analyses are shown in table [T] We see that the two 
Maximum Likelihood methods are in fair agreement in all cases about the significance 
of the correlation, measured by T . However, the number of correlations inferred differs 
strikingly between the two methods for the combined HiRes- AG AS A dataset. Perhaps 
surprisingly, the HiRes method finds more correlations in the AGASA sample alone 
than in the combined dataset. This is a hazard of using the HiRes method that we 
have probed by simulations. Consider the case of combining two sets of events, one 
with poor event resolutions with a number of events aligned to the set of sources, 
and the other set with good event resolutions but all events uncorrelated to the set 
of sources. In such a case the HiRes method significantly underestimates the number 
of correlated events. Conversely, if the aligned events all are taken from the set with 
good resolutions, the HiRes method overestimates the correct number of correlations. 
This explains the discrepancy between the methods for the combined data in table [1] 
since the AGASA events have poor event resolutions relative to the HiRes events. In 
effect, the HiRes method imposes a sort of internal consistency on the distribution of 
correlations, which could be inappropriate, if one data set had systematically incorrect 
a values for instance. 

4.2. BLLacs 



The binned analysis performed in lGorbunov et al.l (J2004 ) on t he sample of 156 BLLacs 



with optical magnitude T^ < 18 from the Veron 10th Catalog (jVeron-Cettv and Veron 



200l[ ) and the 271 HiRes events with E > 10^^ eV showed a correlation at the 10"^ 



level. Applying the HiRes method on this sample we find rig = 8.5 correlations, and 
the fraction T of Monte Carlo runs with greater likelihood than the real data to be 
2 X IO^tJ. Analyzing the same UHECR-BLLac data with the generalized ML method 

J This Us differs from the value (n^ = 8.0) quoted in lAbbasi et al.l ||2006| ) probably due to the HiRes 
exposure map provided to us by HiRes (S. Westerhoff, private communication) being slightly different 
than the one used in the published HiRes analysis. For purposes of comparison of the two methods 
we do not require a perfect exposure map. 
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Table 2. Correlations between BLLacs and UHECRs. The integer 

implementation has been used for the generalized method. 



Generalized ML Method HiRes Method 

UHECR N n(0) n(2) JF uhr :Fhr 

HiRes 271 15 9.2 610""' 8.5 210"" 

AGASA 57 13 6.4 0.04 2.6 0.1 

Both 328 21 12.0 3 10"^ 10.8 110"'' 



we find the number of correlations is n"' — 9.2 with the integer implementation and 
9.3 with the non-integer implementation, and JF = 6 x 10^^ for both. In fact from 
figure [H we see that the HiRes method underestimates ris by 0.7 on average under 
these conditions, so on average this is the expected discrepancy. However as seen from 
figure m the difference between the estimated number of correlations is also consistent 
with the variance from one realization to the next. With la error bars, our best 
estimate is 9.2 ± 2.5 correlations. 

Now including AGASA events too, in Table[2]we list the results for the generalized 
method (integer implementation) and the HiRes method when applied to HiRes and 
AGASA data, both separately and combined. We note that while the calculated 
number of correlations differ between the methods, the difference is within the 
expected dispersion. 

One of the advertised benefits of the generalized ML method is that it provides a 
way to rank sources according to the likelihood of them being correlated with cosmic 
ray events. To exemplify this we list in Table [3] the individual BLLacs and correlated 
HiRes cosmic rays ranked acc ording to the non-integer number of correlations found 
per source (cf. Table 1 in Gorbunov et al.l |2004)). The list includes all cosmic rays 
with positive InTZ, which brings the number of cosmic rays listed to 16. Note that 
the sixth listed cosmic ray is not deemed a correlation by the generalized method, but 
is kept in the list to clarify why the corresponding source is found to have rij > 1 
by the non-integer implementation of the method. Finally, we note that it would be 
desirable to impose restrictions on the sample of BLLacs to include only those within 
an estimated GZK maximum distance for each cosmic ray event, but unfortunately 
the redshift is known for only a fraction of the BLLacs. 

5. Summary and Conclusions 

We have introduced a generalization of the HiRes Maximum Likelihood method, which 
allows the most likely sources of individual events to be identified and ranked. Using 
simulations we have tested the two Maximum Likelihood methods and find that they 
complement each other well: the HiRes method allows a fast way to estimate the 
number of true correlated events, while the new method gives the quality of correlation 
between individual sources and cosmic rays rather than just the total number of 
correlated events. Furthermore, the new method is less sensitive to the validity of 
the estimated angular resolution, and has less systematic bias when candidate sources 
are clustered (as astrophysical sources such as BLLacs are). Of course, for any given 
data set, mere statistical fluctuations can result in conclusions that would not be 
borne out with a larger data set. We conclude that both methods should be used: if 
they disagree markedly on the total number of correlations the data may have some 
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Table 3. HiRes UHECR cross correlations with BL Lacs using the noninteger 
implementation of the generalized Maximum Likelihood method. 



BLLacs 








Ae* (°) 




UHECRs 




Name 


z 


V 


Uj 


InU 


E (EeV) 


RA(°) 


dec (°) 


SBS 1508+561 




17.3 


1.905 


0.83 


3.18 


21.7 


226.56 


56.61 










0.89 


2.90 


11.0 


228.86 


56.36 


MS 10507+4946 


0.140 


16.9 


1.858 


0.60 


4.20 


64.7 


162.61 


49.22 










1.04 


2.08 


24.1 


165.00 


49.29 


GB 0751+485 




17.1 


1.268 


0.32 


4.33 


12.0 


118.77 


48.08 










1.37 


0.21 


15.7 


120.10 


47.40 


lES 1959+650 


0.047 


12.8 


0.991 


0.12 


4.61 


16.4 


300.28 


65.13 


FIRST J11176+2548 


0.360 


17.92 


0.990 


0.12 


4.60 


30.2 


169.31 


25.89 


Ton 1015 


0.354 


16.50 


0.986 


0.37 


4.24 


13.3 


137.22 


33.54 


RXS J01108-1254 


0.234 


17.9 


0.986 


0.37 


4.24 


10.5 


17.79 - 


-12.56 


RXS J03143+0620 




17.9 


0.980 


0.50 


3.89 


27.1 


48.53 


5.84 


RXS J13598+5911 




17.9 


0.970 


0.62 


3.47 


24.4 


210.02 


59.80 


RGB J1652+403 




17.3 


0.958 


0.70 


3.15 


16.5 


253.62 


39.76 


RXS J08163+5739 




17.34 


0.950 


0.74 


2.99 


20.7 


123.78 


56.93 


OT465 




17.5 


0.849 


0.95 


1.98 


10.5 


265.37 


46.72 


lES 2326+174 


0.213 


16.8 


0.527 


1.12 


1.08 


38.1 


352.40 


18.84 



statistical anomaly and be difficult to interpret. 

Although we have tested both Maximum Likelihood methods' ability to reproduce 
the true number of correlations under many different conditions, it is impossible to test 
for every conceivable statistical property of future data sets to which these methods 
could be applied. Thus we recommend, when applying these methods on new data 
samples, that the methods first be tested by simulations using the candidate source 
positions but randomly generated cosmic ray directions with some fraction of them 
Gaussianly aligned to some randomly chosen sources. By repeating this procedure, 
the variance and possible systematic bias in the extracted number of correlations 
can be determined for that particular detector exposure, source and event statistics, 
resolution, etc. 

We applied the ML methods to find correlat ions between x-ray c l usters and 
AGASA events with E > A x 10^^ eV, as done in IPierpaoli and Farrail (|2005l ) but 
without imposing the GZK horizon. We find « 10 excess correlations, with a chance 
probability of 6 x 10~^. The same analysis for HiRes data with E > 10^^ indicates 3—5 
excess correlations, which is expected by chance about 30% of the time. Imposing the 
restriction th at the correlated x-ray cluste r not be farther than the GZK distance cutoff 
was found in IPierpaoli and Farrail ( 20051) to reduce the number of chance correlations 
without significantly decreasing the number of true correlations. Since in this work 
we did not impose the GZK restric tion, we cannot compare the present analysis with 
that of IPierpaoli and Farran ( 20051 ). 

Finally, we applied the generalized Maximum Likelihood method to check 
previously claimed correlations between UHECRs and BLLacs. We corroborate that 
there is a significant correlation b etween BLLac s of the Veron 10th catalog and HiRes 
cosmic rays with E > 10^^ eV (JAbbasi et al.l . 12006 ) . The generalized ML method 
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yields a slightly greater number (sa 9 instead of w 8) excess correlations and a lower 
chance probability (6 x 10~^ instead of 2 x 10^**). Simulations show that the HiRes 
ML method systematically underestimates the number of correlations by 0.7 when 
the source candidates are clustered like the BLLacs, suggesting that the present 
result is the more reliable. The striking difference between the apparen tly significant 
correlation observed in the Northern sky (HiRes) and Auger's null result ( Hararil2007l l 



is a mystery. It may be due to greater average magnetic deflection in the Southern sky, 
or greater incompleteness in the Southern BLLac catalog, or a statistical fluctuation 
or a combination of the above. 
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